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Abstract 

The half lives are calculated for the process of (3^ decay and electron capture for nuclei 
in mass range ~ 65 - 100 relevant for the core of a massive star at late burning stage 
of steller evolution that leads to supernova explosion. The rates of electron capture 
process are also computed for nuclei relevant during presupernova stage of massive 
stars. These half lives and rates are calculated by expressing [3^ Gamow- Teller decay 
strengths in terms of smoothed bivariate strength densities. These strength densities 
are constructed in the framework of spectral averaging theory for two body nuclear 
Hamiltonian in a large nuclear shell model space. The method has a natural extension 
to weak interaction rates for r and rp-processes. 
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1 Introduction 

The late evolution stages of massive stars (> 8M Q ) is strongly influenced by 
the weak interaction processes. In the pre-collapse stage of a star the important 
weak processes are electron captures and beta decays. These weak processes 
influence the value of Y e , the electron fraction. As Chandrasekhar mass is 
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proportional to Y^, the value of Y e is crucial for the collapse of a supernova star 
and the subsequent evolution. The stability (or collapse) of the star at late stage 
is determined by among other things, the balancing of the gravitational collapse 
by electron degenerate pressure. While electron capture reduces the number 
of electrons for pressure support, f3~ decays increases the same. Therefore 
these rates are important ingredients for building a model for the core collapse 
supernova. With this in view, in this work we have investigated the rates 
and half lives for weak interaction processes for some of the nuclei relevant for 
presupernova and supernova environment with 65 < A < 100. We first calculate 
the weak interaction strength densities using a theory based on principles of 
statistical nuclear spectroscopy - spectral averaging theory. 

Among the earlier calculations for weak interaction rates (/3 ± , EC etc.), 
Fuller, Fowler, Newman |IJ for the first time dealt with these rates at finite 
temperature in the mass range 21 < A < 60. Aufderhide et al |3| emphasised 
the need for calculating weak interaction rates for A > 60 nuclei as presupernova 
stars generate neutron rich cores with mass A > 60. Aufderhide et al and later 
Kar et al [f| have made rate calculations in the mass A > 60. In this work, 
we have explicitly constructed the bivariate weak interaction strength densities 
for calculation of weak interaction half lives and rates. As in an earlier work, 
using these methods, f3~ decay rates have been calculated || for nuclei in the 
regime 55 < A < 65 and extended to EC for a few fp-shell nuclei, in this work 
we have computed EC half lives for the nuclear mass region 65-81 that are 
still relevant for supernova stars using smoothed form for EC strength densities 
in the framework of spectral averaging theory. We also calculate rates for some 
sample nuclei in this region and this would be extended to neutron rich nuclei 
later. 

After the spontaneous nucleosynthesis ceases in late stage of steller evolu- 
tion, heavier elements in the star are synthesised by r-process, s-process and 
rp processes. The drip line nuclei play an important role then for the under- 
standing of stellar processes after the star runs out of fuel. With this in view 
we have calculated the /3 ± half lives for the drip line neutron nuclei in the mass 
region 81 < A < 100. In this case too the calculations are made by constructing 



2 



bivariate strength densities for /3 ± operator under the framework of spectral 
averaging theory. 

The spectral averaging theory which is based on statistical nuclear spec- 
troscopy, was started with Bethe's statistical mechanical level density formula, 
Wigner's treatment of spectral fluctuations using matrix ensembles and French's 
embedded ensembles and Gaussian densities. The smoothed forms of spectro- 
scopic observables follow from the action of Central Limit Theorems (CLT) in 
nuclear shell model spaces. The statistical spectroscopy is in deriving and ap- 
plying the smoothed forms in indefinitely large spaces with interactions by using 
unitary group decompositions (of Hamiltonians and the spectroscopic spaces), 
CLT's locally, and convolutions - the resulting theory is the Spectral Averaging 
Theory in Large Shell Model Spaces (SAT-LSS) ||, [7|, ||, ||, [10|, |TTjl - Here it is seen 
that the essential role of interactions is to produce local spreadings of the non 
interacting particle (NIP) densities and the spreadings are in general Gaussian 
in nature. The spectral averaging theory in large shell model spaces has im- 
portant nuclear physics applications like calculations of nuclear state and level 
densities, occupation and spin cut-off densities and calculations of occupation 
numbers and spin cut-off factors fTU|| . 

The theory has been extended to calculate the smoothed form for interacting 
particle (IP) bivariate transition strength densities (I^ (Ei,Ef)) for a one plus 
two body nuclear Hamiltonian H and a transition operator O. This bivariate 
strength density (Iq (Ei, Ef)) takes a convolution form [^, [12], [13] with the non 



interacting particle (NIP) strength densities being convoluted with a spreading 
bivariate Gaussian due to irreducible two body part of the interaction. Math- 
ematically, 1q =>- 1q <S> Po, where I h is the strength density due to effective 
one body part h of the interacting Hamiltonian and is a zero centred bi- 
variate Gaussian due to the irreducible two body (off-diagonal) part V of the 
interaction. In spectral averaging theory the NIP part 1^ are constructed by 
calculating a few lower order moments of 1^ and explicit analytical formulae for 



NIP strength densities are worked out in ||13| . This method termed as moment 
method is a convenient way for rapid construction of NIP strength densities. 
The IP bivariate strength densities thus constructed has been used to calculate 
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(3 decay rates for some fp-shell nuclei around the mass range 55 < A < 65 



Giant Dipole Resonance (GDR) cross-sections |14| etc. 

The smoothed form for IP strength densities for (3^ decay and EC Gamow- 
Teller transition operator can therefore be constructed using the formalism of 
spectral avearaging theory in a large shell model space. With these, the half 
lives and rates for electron capture processes at finite temperature relevant for 
presupernova environment and half lives for (3^ decay processes are calculated 
for the nuclei mentioned earlier. 

The paper is organised as follows. In Section 2 we give the formalism. 
Section 3 deals with the calculational methods that includes choice of shell 
model space, the single particle energies (s.p.e.) etc. The results are discussed 
in Section 4. In section 5 we conclude with some discussions and future outlook. 



2 Formalism 

2.1 Unitary decompositon of nuclear Hamiltonian 

Given m number of particles distributed in a shell model space we have an m 
- particle shell model space. Each shell model orbit or single particle orbit a 
with degeneracy N a = 2j a + 1 is called a spherical orbit. A group of spherical 
orbits is called a unitary orbit OL . In what follows we will use a, (3 etc. (smaller 
in size) to denote spherical orbits and OL , (3 etc. (larger in size) for unitary 
orbits. Thus we have for m number of particles distributed in this shell model 
space, spherical configuration m = m a ,mp^... (m a , mp etc. are number of 
particles in spherical orbits a, (3 etc. respectively) and unitary configuration 
[m] = m a ,mg , .... Using the convention that the spherical orbits a belongs 
to unitary orbit OL the number of single particle orbits in uniary orbit OL is 

A further decomposition of m-particle space is possible by attaching an s a 
label to each spherical orbit a where s a for lighter nuclei denotes sTiuo excitation 
value. With this the m- particle space can be decomposed into S'-subspaces as 
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follows: 

m^E^^^EWH^E^^E^ (!) 

One can recognise here the appearence of U(N) group which is generated 
by the N 2 operators a] m aj ama ; a, (3 = 1,2, ...,N. For m identical particles 
I N \ 

therefore there are antisymmetric states forming an irreducible repre- 

V m / 

sentation (irrep) for the group C/(iV) usually denoted by Young shape {l m }. 
The only scalar operator is the number operator n as it remains invariant under 
the transformation of the U(N) group. A given operator can be decomposed 
into tensor operators (belonging to a definite irrep of U(N)) with respect to 
U(N) group. 

Thus the unitary group U(N a ) acting on each spherical orbit a generates 
m a of spherical configurations m, i.e. m behaves as {l m «} © {l" 1 / 3 } (g) .... with 
respect to the direct sum group U(N a ) © U(Np) © ... and the scalar operators 
are m Q 's. Therefore for a given nuclear two body Hamiltonian H = h(l) + V(2), 
it should be obvious that non interacting particle part h(l) of if is a scalar with 
respect to spherical configuration group and 

h = J2^n a = h^ (2) 

where e a is the single particle energy (s.p.e.) of the spherical orbit a. For the 
spherical configuration scalar part (must be a second order polynomial in 
n a ) of residual interaction V(2), one has 

n a {np - 5 a p) 

In the above V a p is the average two body interaction given by 



( 1 + Oa/3) 



V aP = {N^}- 1 ^ J + 1)^(1 + 5 aP )} (4) 
j 

where N a p = N a (Np — 8 a p) and J is the angular momentum. The remaining 
non-diagonal part V = V(2) — V' ' of V(2) is an irreducible two body part as 
there cannot be an effective one body part of V(2) with respect to spherical 
configuration group. 
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The same idea can be applied for the decomposition of Hamiltonian under 
unitary configuration group. In this case, the unitary group U(N a ) acting 
in each unitary orbit Ct generates m« of unitary configuration [m]; i.e. [m] 

m n 

behaves as {l m a }©{1 P }©.... with respect to the direct sum group U(N a )ffi 
U(Np )©... and the scalar operators are ma 's. The decomposition relevant for 
the present is the spherical scalar part + of the Hamiltonian. They will 
have [0] © [1] © [2] tensorial parts with respect to unitary configuration direct 
sum group. 

Thus after the unitary decomposition of the nuclear Hamiltonian H we have 
H = h(l) + V(2) 

^ fc[0][0] + ^[0][l] + y[0][O] + y[0][l) + y[0][2] + y ^ 

This has been demonstrated that the contribution of V^ 2 ! part is small 
(< 5%) all across the periodic table by calculating its norm for the case of 
ds, fp, 10-orbit and 15-orbit interactions and phenomenological interactions like 
surface delta interaction and pairing + Q.Q interactions [0, |HJ. Therefore I/MI 2 ! 
can be negleced for all practical purposes. 

Thus Eq. 5 reduces to 

H = h m + h m] +V m + V [m +V 

= h + V (6) 

where h is the effective one body part of H and V is irreducible two body part. 

For a unitary configuration [m], the unitary decompositions of nuclear Hamil- 
tonian H are given by 



a \aea / 



h [0][l 



y[o][o] 



a >p 1 + « P 1 



V, 



a P 



N 



a P 
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yIPJP 




-<5 



> n n 



E(^ 



a/3 



S a p)V a p \x{N p -25 a/3 }- 1 (7) 



In the above and for rest of the calculations we consider S'-conserving part (see 
discussions above Eq. 2 for S quantum number) of the interaction. The S- 
mixing part Vs-mix of ^(2) represents admixing between distant configurations 
(at least 2Tiuj away from each other) and this leads to multimodal form of densi- 
ties ]n| unlike the unimodal forms. Moreover GT/? 1 * 1 operator does not connect 
different S-subspaces. Hence the Vs-mix is not considered for the rest of the 
calculations. 

Given H = h + V the IP strength density lQ =h+v (Ei, Ej) for a transition 
operator O will take a bivariate convolution form ]12[ with the two convoluting 
functions being NIP strength densities Iq and a normalised spreading bivariate 
Gaussian Po-biv-g due t° V (interactions), 



T H=h+V 

1 o 



1 o 



V 

PO;BIV-G 



[E h Ej 



(8) 



For our case, O = 0(GT). In large spectroscopic spaces with protons and 
neutrons (pn), this GT bivariate strength density can be partitioned in different 
unitary configuration subspces and S'-subspaces (Eq. 1) and can be written as 
(identifying that 0(GT) does not connect two different subspaces) 



l O(GT) 



In geneal for a transition operator O and for a nuclear Hamiltonian if, the 
strength density \ H '>( m > m ) [ s given as 



E 

s [m*,mk],[m£,m£]eS 



lH ;(m,m') = i™\ E ')\(E' m'\0\E M)\ 2 I m (E) 
= ((0^6{H - E')OS{H - E))) m 



(10) 
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where I m ' (E') and I m (E) are final and initial state densities and (( )) m repre- 
sents a trace over m particle space. By the action of CLT in the spectroscopic 
space, the density will be a bivariate Gaussian and it is demonstrated in |10|. |l| 
by constructing the exact NIP strength densities (and its s-decomposition) for 
0(GT) operator and then comparing this with the smoothed Gaussian form. 

The smoothed form for strength density is constructed with marginal cen- 
troids 6\, €2 and variances o\ , a\ and calculating few lower order central bivariate 
moments given by 

e 1 = (OWH) m /(0^0) m 

e 2 = (0 ] HO) m /(0 ] 0) m 

a\ = (O^OH 2 ) m /(0 ] 0) m 

a 2 = (0 ] H 2 0) m /(0 ] 0) m 




Therefore the calculation of moments are in fact calculation of m-particle av- 
erages or traces of the operator in question. This is done by first calculating 
a few basic traces and then propagating them in m-particle space. As here we 
are dealing with unitary configurations [m] and unitary configuration densities 
(Eq. 9), we would require the unitary configuration traces of the type ( )[ m '. 
These moments M pg ([m]) for the construction of NIP strength densities Vq for 
a one body transition operator O with p + q < 2 are calculated in details in 
Ref. |L0|, [|. As the GT operator is of the type 0(GT) = e a pa) a ap (e a p is the 
single particle matrix elements), for a given initial configuration [nij], the final 
configuration [mj] then is obtained uniquely as [my] = [m^] x (l^l^j. Thus one 
also obtains the partial moments M M ([mj], [mj]) for I h . For a proton neutron 
(pn) configuration one writes [m] as [m p ,m n ]. 

For the construction Gaussian spreadings p ' p ' n ' p ' n (GT)(x,y), in 
Eq. 9, the following approximations are adopted. The marginal centroids, 
M 10 = (O^OV) m p> m ™/(0^0) m p' m ™ ~ (V) m p' m « = 0, as V is traceless; M 01 = 

f f f f f f 

(O^V~0) mp,mn / (O^ 0) mp ' mn ~ (v) m P' mn = and the marginal variances given 
by traces of the type (O^OV 2 ) and (OW 2 0) are equal to (V 2 ). 
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2.2 Formalism for ^ decay half lives and electron cap- 
ture rates at finite temperature 

The weak interaction {(3^ and electron capture (EC) for the present case) 
rate T(Ei — > Ef) is the number of weak processes per second from a given 
initial state \Ej) of the parent nucleus to the final nuclear state \Ef) and 
T(Et -> E f ) oc [gyB F (Ei -> Ef) + g\B GT (Ei -> £?/)], where gv and g A are 
respectively the vector and axial vector coupling constants and Bp and Bqt are 
the Fermi and Gamow - Teller transition strengths respectively. Including the 
phase space factor / that incorporates the dependence of the rate on nuclear 
charge Ze and the available energy for the weak process under consideration, 
T takes the form T{E, t -> E f ) = Cf[g^B F (E t -> E f ) + g\B GT (E t -> 
where C is a constant with Q the GS Q- value of the weak process ( / 9 ± or EC), 
Qi = Q + Ei. The value of C is fixed by using the log ft values from pure 
Fermi transitions. One can write down the expressions for ground state half 
lives and f3 ± decay and EC rates at finite temperature. For present beta decay 
calculations we have neglected Fermi term Bp as the Fermi strength is concen- 
trated in narrow domain and high up in energy (the centroid ~ 1AAZ A^ 1 ^ 3 
MeV, width ~ 0.157Zy4 -1 / 3 MeV). For (3 + decay and electron capture pro- 
cesses, Fermi transitions are not possible for nuclei with N > Z, since for 
mother isospin T (= (N — Z)/2), the isospins of the states of daughter are 
greater by 1, i.e. T daughter = \{(N + 1) - (Z - 1)}/2| = T + 1. Writing 
Ba^JiEi -> JfEf) = (2J l +l)- 1 E\(EfJfMf\(0 GT )';\E l J l M l ) 2 which in contin- 
uous version becomes {/(-^i)}" 1 Eaes^eE/^ K-E/a|(0ar)J|-Ei/3)| 2 , the expres- 
sion for GS half life is 



' 9 -A 2 3£ 
,9v 



ti/ 2 {GS) = {6250(s)}x 

l O(GT)( E GS,Ef) 



I H (E 



GS 



f(Z)dE f (12) 



In the above equation £ is the so called quenching factor and the factor 3 comes 
because of the definition of l h . The usual values of £ is 0.6 for f3~ decay and 

2 



0.5 for EC//3 + decay and ( — -J = 1.4 as given in [1?|. The weak interaction 
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rates are given by 
A(T) = hl2iS " 



6250 



J dE i e~ Ei < /kBT I H (E i ) J^dE 



X 

\9v. 



3£ } x 



I H {E„ 

£n2(s- r 
6250 



/(Z,T) 



e- E ^ BT I H {E i )dE i 



i -l 



jdE t e- E *' kBT J Qi dE f 



9a 
.9v. 



3£ )■ Iq{gt) 



(E u E f )f(Z,T) 



(13) 



The phase space integral / in the above eqautions are given by, for the case of 
EC @ 

/ = / dEE(E 2 -m 2 e ) 1 / 2 (E-Q l ) 2 F c (Z,E)x 



1 + exp[(.B - n e )/k B T] 



1 



1 + exp[(i£ — /x n «) / fc B T] 



(14) 



where E m i n is greater among —Qi and m e c 2 . In the actual calculation the 

integrand is expressed in terms of electron mass m e and the integration limits 

are also accordingly changed. In the above k B stands for Boltzmann constant, 

He is the chemical potential for electron and fi u is the same for neutrino. For 

the case of (3~ decay the phase space factor reads as 

-oF c (Z,6 e )6 e (e2-l)V2 (eo _ e ^ 2 



/ 



-de P 



(15) 



h {1 + exp[(/j, e - e e )/(k B T) e }} 
In the above equation e = E /m e ; (E = Qi — Ef), /i e = /i/m e and {k B T) e = 
k B T/m e . F C (Z,E) in the above equations give the Coulomb factor and for EC 
is given by |18[ 

F C (Z,E) = F'{Z)F' e {E) 

87raZ(l + S) (2R\ 2S - 2 < 1 + 0.577(aZ) 2 \ 



Fc |r(l + 2s)| 2 V h J 

F' e (E) = (E 2 - m 2 ^- 1 ' 



2 + 0.577(aZy 



E 



(E 2 - m 2 ) 1 / 2 



1 - e 



-2-Kb 



(16) 
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In the above, a is the fine structure constant. The value of S in the present 

i 



l-e- 2 *" 6 



calculation is assumed to be 1 (T^j and following [JTSj the term 
expressed as (27i9) and for Qi < m e c 2 

{2,9) = 2naZ eX ^- 1 / / ^ (l + 2\ + 2\>) 
X = kBT/m e c 2 



where K 2 is second order modified Bessel function. For Qi > m e c 2 we use 



is 



;i7) 



(2710) 



2naZ- 



m 



2)1/2 



(18) 



For f3 process the Coulomb factor is taken as given in Schenter and Vogel 
| and the expression reads as 

F c (Z,e) 6 
a(Z) 



exp[a(Z) + /3(Z)Ve 



HZ) 



-0.811 + 4.46(-2)Z + 1.08(-4)Z 2 
-8.46(-2) + 2.48(-2)Z + 2.37(-4)Z 2 
0.673 - 1.82(-2)Z + 6.38(-5)Z 2 
1.15(-2) + 3.58(-4)Z - 6.17(-5)Z 2 



(e- 1) < 1.2 
(e- 1) > 1.2 
(e- 1) < 1.2 
(e - 1) > 1.2 



(19) 



(20) 



For the electron chemical potential we have used theexpression given in 

In the above p-j is the matter density p in units of 10 gms/cc, T is the tempera- 
ture T expressed in MeV and Y e is the electron fraction. We put p u , the neutrino 
chemical potential to be zero as for the densities we consider the neutrinos are 
free streaming. 



3 Calculational Procedure 

For the construction of strength densities we have selected a 9-orbit shell model 
space both for proton and neutron with 56 Ni as core, consisting of the spherical 
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orbits 3p 3/2 , 3/ 5 / 2 , 3pi/ 2 , 4gr 9 / 2 , 4d 5 / 2 , 4# 7/2 , 4si/ 2 , 4d 3/2 and 5/tn/ 2 . The s.p.e. in 
MeV are 26.58, 26.19, 29.09, 33.91, 38.52, 42.47, 42.30, 43.15, 58.44 respectively. 



The initial values of s.p.e's are taken from ref |f20[ . The renormalisation effects 
due to the closed core (in this case s shell, p shell, ds shell and f 7 / 2 orbit) are 
then incorporated. This effect not only renormalises fp — sdg and sdg — h\i/2 
separation but also renormalises the single particle energies. They can be eval- 
uated from V a/ 3 and V a g discussed in Sect. 2.1. The unitary orbits are {3p 3 / 2 , 
3/ 5 /2, 3p 1/2 }, {4sf 9/2 }, {4c/ 5/2 , 45( 7/2 , 4s 1/2 ,4d 3/2 }, {5/i 11/2 }. Thus each of the pro- 
ton and neutron shell model space has been divided into four unitary orbits. For 
the 2-body residual interaction, we have used a phenomenological interaction 
namely pairing + Q.Q interaction [21] with the strength \ = 242/A -5 / 3 . 

The nuclei chosen for different rate and half life calculations are given in 



Table 1. These nuclei have been chosen from Ref. |22] 

Using the formalism given in Sect. 2 one can now construct the weak in- 
teraction strength densities (bivariate Gaussian form). For the state densities 
I H (E), required for calculation of half lives or rates (Eq. 12, 13), we adopted 



the formula given in Dilg et al |p3]| . This has been demonstrated in [|H], |11 
that state densities I H {E) obtained from spectral averaging theory represent 
very well the results obtained from Dilg et al state density formula. Besides the 
assumption of marginal centroids and variances (Sect. 2.1), we also assume the 
bivariate correlation coefficient £ to be independent of configuration. Moreover, 
for two-body EGOE (Embedded Gaussian Orthogonal Ensemble) ||, one has 
the result ( ~ 1 — 2/m, where m is the number of active nucleons. Therefore, 
in the present calculation, ( is first evaluated by treating it as a parameter with 
the EGOE form for ( = a + b/m and a and b is evaluated by calculating f3~ 
decay half lives for 18 nuclei in the mass range considered and then comparing 
them with the experimental values of the same. 

To this end, we fix a and b by calculating £ for various values of a and b 
with the constraint that the value of ( lies within 0.6 and 0.9. The value of 
bivariate variance (approximated as state density variance, (V 2 ) (Sect. 2.1)) is 
taken to be 16.5 MeV 2 . The values of a and b are found out by minimizing the 
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quantity ^ {£og(T[/ 2 ) ca i — £og(T{/ 2 )expt) 2 - The values of a and b are obtained 

i=nuclei 

respectively as 0.77 and 0.3. The values of ( and j3~ half lives for the nuclei 
considered are listed in the Table 2. 

Having obtained the values of the parameters a and b for calculation of 
we now proceed to calculate f3 + decay half lives for the nuclei given in Table 1. 
For this calculation, we keep ( unchanged from the value obtained for f3~ decay 
half life calculations. But here we have parametrised the variance as = A + 
B/m and obtained the parameters A and B by minimising {^°g{ T i/2)cai ~ 

i=nuclei 

£og(rl/ 2 ) expt ) 2 for the four nuclei considered for (3 + decay (Table 1) much the 
same way as determination of ( in f3~ decay calculations. For calculating [3 + 
half lives we use Eqs. 12 and Eqs 16 - 18. The values of A and B are obtained 
as A = 20.5 and B = 1.9. The values of calculated half lives and the variance 
o"y for the four nuclei considered for (3 + decay are shown in Table 3. 

Calculation of electron capture (EC) half lives are similar to that of f3 + 
calculations. We use the parameters a and b obtained from (3~ calculations 
to obtain correlation coefficients ( for the five nuclei (Table 1) chosen for EC 
calculations. The variance ay is parametrized as Uy = A + B/m and the values 
of A, B and consequently cry are found out by minimising (£og(Ty 2 ) ca i — 

i=nuclei 

£og(T{/ 2 ) expt ) 2 . In this case the values of A and B are obtained as A = 11.5, 
B = —4.6. For these calculations we have used the Eqs. 12, 14, 16 - 18. The 
calculated half lives and are shown in Table 4. 

We have also calculated EC decay rates for typical densities and tempera- 
tures relevant for presupernova stars. The rates are evaluated for three different 
temperatures (3 x 10 9 , 4 x 10 9 and 5 x 10 9 in °K), three different stellar densities 
(10 9 , 10 8 and 10 7 gms/cc) and for three values for electron fraction Y e (0.50, 
0.47, 0.43). The calculations use Eqs. 13,14, 15-18 and Eq. 20 are used. The 
calculated EC rates alongwith temperatures, Y e and stellar densities for each of 
the five nuclei considered, are given in Table 5. 
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4 Discussions and Conclusions 



We have calculated the half lives and rates for weak interactions for a number 
of nuclei spanning a wide range (65 < A < 100) in the periodic table. These 
calculations are performed by explicitly constructing smoothed form for bivari- 
ate strength densities for the weak interaction operator (in this case Gamow 
Teller operator) using convolution form within the framework of spectral av- 
eraging theory in nuclear physics. The lower mass region of the above range 
is relevant for presupernova and supernova stars and the higer mass region is 
relevant for r-process nucleosynthesis that synthesises heavier nuclei. We have 
used here, the principles of spectral averaging theory for calculation of EC rates 
for different temperatures, densities and electron fraction values in presupernova 
environment. Earlier such calculations with spectral averaging theory have been 
performed for (3~ decay rates for presupernova stars and for a few fp-shell nuclei 
pijj . The (3 + and (3~ decay half lives are calculated for the nuclei in the range 



A > 80. The calculation of the rates for typical r-process conditions for these 
processes will be taken up in near future. 



The authors thank A. Ray for helpful discussions. 
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Table Captions 



Table 1 Nuclei chosen for different weak processes. 

Table 2 Calculated and experimental f3~ decay half lives and correlation coef- 
ficient (. Method of calculating ( is given in the text. 

Table 3 Calculated and experimental f3 + decay half lives and variance a^. 
Q-values are also given. 

Table 4 Calculated and experimental EC half lives and variance ay. Q-values 
are also given. 

Table 5 Electron Capture rates for 65 Ge, 69 Se, 73 Kr, 77 Sr, 81 Zr. 
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Table 1 



Weak Interaction 


Nuclei Chosen 


Electron Capture 


65 /nr_ 69 q„ 73 ^„ 
32 ( - Te 33> 34 35 ? 36 Ar 37, 

3s*S'r39,?i^r4i 


(3 + decay 


42MO43, 44-RM45, 46-Pc?47; 

HCd A9 


(3~ decay 


30^51) 31^50, 33^548 
32 (jre 50, 31^052, 33^50 
32^ e 52, 34 Se 5 \, 33AS54 
37-R65I, 345*655, 35-8^54 

36 iir 54, 37-KO54, 35-D^57 

92 a„ 93 94 v 
38°' 54, 36-"- ' 57, 39*55 
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Table 2 



Sr No. 


Nucl. 


z 


Qval 

(MeV) 


rpexpt. 
1 1/2 


rpcalc 
1 l/2 

/ 


c 


1 


81 


30 


11.9 


0.29 


0.819 


0.7820 


2 


81 Ga 


31 


8.32 


1.221 


1.483 


0.7820 


3 


81 As 


33 


3.86 


33.3 


3.662 


0.7820 


4 


82 Ge 


32 


4.7 


4.6 


53.088 


0.7815 


5 


83 Ga 


31 


11.5 


0.31 


0.127 


0.7811 


6 


83 As 


33 


5.46 


13.4 


21.134 


0.7811 


7 


84 Ge 


32 


7.7 


1.2 


2.415 


0.7807 


8 


85 Se 


34 


6.182 


31.7 


106.261 


0.7803 


9 


8 Ms 


33 


10.3 


0.73 


0.0632 


0.7797 


10 


88 M> 


37 


5.13 


1066.8 


27227.014 


0.7794 


11 


89 Se 


34 


9.0 


0.41 


0.117 


0.7791 


12 


89 £> r 


35 


8.16 


4.4 


1.465 


0.7791 


13 


90 


36 


4.39 


32.32 


413.661 


0.7788 


14 


91 


37 


5.86 


58.4 


290.453 


0.7786 


15 


92 5r 


35 


12.2 


0.343 


0.096 


0.7783 


16 


92 Sr 


38 


1.67 


9756. 


110696.591 


0.7783 


17 




36 


8.6 


1.286 


3.706 


0.7781 


18 


94 y 


39 


4.0 


67320. 


2590.043 


0.7779 
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Table 3 



Nucleus 


z 


Qp+ 

(MeV) 


rpexp 
1 l/2 

sec 


r-pCalc 
1 l/2 

sec 


<4 
MeV 2 


85 Mo 


42 


8.6 


0.4 


0.076 


20.566 


89 Ru 


44 


8.92 


0.3 


0.22 


20.558 


93p d 


46 


9.5 


0.3 


0.14 


20.551 


97 Cd 


48 


9.6 


0.2 


0.32 


20.546 



Table 4 



Nucleus 


Z 


Qec 
(MeV) 


rpexp 
1 l/2 

sec 


rpCalc 
J l/2 

sec 


a 2 
MeV 2 


65 Ge 


32. 


6.24 


30.9 


12.79 


10.989 


69 Se 


34. 


6.78 


27.4 


10.32 


11.146 


73 Kr 


36. 


6.65 


27. 


12.58 


11.229 


77 Sr 


38. 


6.85 


9. 


9.35 


11.281 


81 Zr 


40. 


7.16 


15. 


14.04 


11.316 
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Table 5 



Nucleus 


p(gms/cc) 






Temperature 


in 


°K 








Ye 


3 x 10 9 


4 x 10 9 


5 x 10 9 












Rates (s 


V ) 






65 Ge 


10 9 


0.50 


2.77 x 10" 


-3 


2.89 x 10" 


-3 


3.07 x 10" 


-3 






0.47 


2.56 x 10" 


-3 


2.68 x 10 


-3 


2.86 x 10" 


-3 






0.43 


2.28 x 10" 


-3 


2.41 x 10" 


-3 


2.58 x 10" 


-3 




10 8 


0.50 


2.05 x 10" 


-4 


2.47 x 10" 


-4 


3.13 x 10" 


-4 






0.47 


1.92 x 10" 


-4 


2.34 x 10" 


-4 


2.98 x 10" 


-4 






0.43 


1.76 x 10" 


-4 


2.16 x 10" 


-4 


2.79 x 10" 


-4 




10 7 


0.50 


2.96 x 10" 


-5 


5.67 x 10" 


-5 


1.10 x 10" 


-4 






0.47 


2.86 x 10" 


-5 


5.56 x 10" 


-5 


1.10 x 10" 


-4 






0.43 


2.71 x 10" 


-5 


5.43 x 10 


-5 


1.09 x 10" 


-4 


69 Se 


10 9 


0.50 


2.81 x 10" 


-3 


2.92 x 10 


-3 


3.10 x 10" 


-3 






0.47 


2.60 x 10" 


-3 


2.72 x 10 


-3 


2.89 x 10" 


-3 






0.43 


2.32 x 10" 


-3 


2.45 x 10" 


-3 


2.62 x 10" 


-3 




10 8 


0.50 


2.20 x 10" 


-4 


2.64 x 10" 


-4 


3.31 x 10" 


-4 






0.47 


2.07 x 10" 


-4 


2.50 x 10" 


-4 


3.16 x 10" 


-4 






0.43 


1.89 x 10" 


-4 


2.31 x 10 


-4 


2.96 x 10" 


-4 




10 7 


0.50 


3.26 x 10" 


-5 


6.16 x 10" 


-5 


1.19 x 10" 


-4 






0.47 


3.14 x 10" 


-5 


6.05 x 10" 


-5 


1.18 x 10" 


-4 






0.43 


2.99 x 10" 


-5 


5.90 x 10 


-5 


1.17 x 10" 


-4 


73 Kr 


10 9 


0.50 


2.18 x 10" 


-3 


2.27 x 10 


-3 


2.40 x 10" 


-3 






0.47 


2.01 x 10" 


-3 


2.11 x 10" 


-3 


2.24 x 10" 


-3 






0.43 


1.80 x 10" 


-3 


1.90 x 10 


-3 


2.03 x 10" 


-3 




10 8 


0.50 


1.72 x 10" 


-4 


2.05 x 10- 


-4 


2.58 x 10" 


-4 






0.47 


1.61 x 10" 


-4 


1.95 x 10" 


-4 


2.46 x 10" 


-4 






0.43 


1.48 x 10" 


-4 


1.80 x 10" 


-4 


2.30 x 10" 


-4 




10 7 


0.50 


2.55 x 10" 


-5 


4.82 x 10" 


-5 


9.24 x 10" 


-5 






0.47 


2.46 x 10" 


-5 


4.73 x 10- 


-5 


9.19 x 10" 


-5 






0.43 


2.34 x 10" 


-5 


4.61 x 10 


-5 


9.13 x 10" 


-5 
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Table 5. Contd. 



Nucleus 


p(gms/cc) 






Temperature 


in 


°K 








Y e 


3 x 10 9 


4 x 10 9 


5 x 10 9 












Rates (s 


r ) 






77 Sr 


10 9 


0.50 


1.86 x 10" 


-3 


1.90 x 10 


-3 


1.98 x 10- 


-3 






0.47 


1.72 x 10" 


-3 


1.77 x 10" 


-3 


1.86 x 10" 


-3 






0.43 


1.54 x 10- 


-3 


1.60 x 10 


-3 


1.68 x 10- 


-3 




10 8 


0.50 


1.56 x 10" 


-4 


1.84 x 10" 


-4 


2.26 x 10- 


-4 






0.47 


1.47 x 10- 


-4 


1.74 x 10 


-4 


2.16 x 10- 


-4 






0.43 


1.35 x 10- 


-4 


1.62 x 10 


-4 


2.03 x 10- 


-4 




10 7 


0.50 


2.39 x 10- 


-5 


4.41 x 10- 


-5 


8.27 x 10- 


-5 






0.47 


2.30 x 10- 


-5 


4.33 x 10 


-5 


8.23 x 10- 


-5 






0.43 


2.19 x 10- 


-5 


4.22 x 10 


-5 


8.18 x 10- 


-5 


81 Zr 


10 9 


0.50 


1.45 x 10- 


-3 


1.48 x 10- 


-3 


1.55 x 10- 


-3 






0.47 


1.34 x 10- 


-3 


1.38 x 10" 


-3 


1.44 x 10- 


-3 






0.43 


1.20 x 10- 


-3 


1.24 x 10- 


-3 


1.31 x 10- 


-3 




10 8 


0.50 


1.15 x 10- 


-4 


1.36 x 10" 


-4 


1.68 x 10- 


-4 






0.47 


1.09 x 10- 


-4 


1.29 x 10 


-4 


1.60 x 10- 


-4 






0.43 


9.96 x 10- 


-5 


1.19 x 10 


-4 


1.50 x 10- 


-4 




10 7 


0.50 


1.72 x 10- 


-5 


3.20 x 10 


-5 


6.04 x 10- 


-5 






0.47 


1.66 x 10- 


-5 


3.14 x 10- 


-5 


6.01 x 10- 


-5 






0.43 


1.58 x 10- 


-5 


3.06 x 10 


-5 


5.97 x 10- 


-5 
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